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ABSTRACT 


A fixed base helicopter simulation was implemented 
which may be easily altered to represent a variety of con- 
ventional single rotor aircraft. Test subjects performed 
90-second tracking tasks utilizing two cockpit displays for 
longitudinal control during landing approach. The performance 
data generated were compared with predicted values from an 
optimal pilot modeling technique. Tentative conclusions 
were drawn as to the effectiveness of each display as well 


as the predictive power of the pilot model. 
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LIST OF SYMBOLS 


g Acceleration of gravity, ft/sec 
ii Aircraft mass moment of inertia about y stability axis, 
y 2 

Slug-ft, 


M Moment about y stability axis, ft-lbs. 


m Aircraft mass, slugs. 
q Aircraft angular rate about y stability axis, radians/sec. 
Un Aircraft reference flight velocity, ft/sec, 
u Perturbation flight velocity along x stability axis, ft/sec. 
oe Horizontal wind gust velocity, ft/sec. 
Ww Perturbation flight velocity along z stability axis, 
ft/sec. 
We Vertical wind gust velocity, ft/sec. 
X Force component along x stability axis, lbs. 
Z Force component along z stability axis, lbs. 


- Cyclic pitch control input, ft. displacement at the 
pilot's hand. 


6 Collective pitch control input, ft. displacement at 
tne pilet*s hand. 


C) Perturbation pitch angle about y axis, radians. 


The following stability derivatives are defined for a 


reference flight condition in stability axes. 
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I, INTRODUCTION 


The use of vertical take-off and landing aircraft has 
increased significantly in both military and civil applica- 
tion. Helicopter flight under instrument flight rules has 
been accomplished on a mOneine basis for quite a while but 
low speed and terminal operations still require visual 
references outside the aircraft. The military services have 
had some success in low speed instrument flight but it is 
generally conceded that much work needs to be done in the 
area , information SSUES One popular approach 
has Beenie lise oe aucetean i e tio lave utilizing cathode- 
ray-tubes which have the advantage of almost unlimited 
flexibility in format and symbolism (References 1 and 2). 
Additionally, extensive work has been done in the field of 
pilot modeling and pilot-in-the-loop system performance (Ref.3); 
but in order to optimize pilot performance it is necessary 
to consider the nature and extent of information presented 
to the pilot. \ this project was aimed at the construction 
of a fixed base helicopter flight simulator and the comparison 


of real pilot performance with pilot model predictions under 


a limited range of conditions and using two visual displays. 
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II, METHOD OF INVESTIGATION 


A helicopter simulation was constructed to represent 
the longitudinal flight dynamics of an unaugmented (i.e., 
without the usual mechanical stabilization bar) UH-1H 
helicopter. A fore and aft moving control stick provided 
cyclic pitch control inputs and a collective pitch level 
produced power commands. A hybrid computer gave real time 
solutions to the helicopter equations of motion, generated 
each of two cockpit displays, and computed performance 
data. Performance values for each display mode were 
analytically predicted using optimal pilot modeling tech- 
niques and compared with simulation data for two test 


subjects. 


III, APPARATUS AND PROCEDURES 


The longitudinal motion of a helicopter can be described 


by the following equations of motion: 
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Sy represents cyclic control input, os collective input, 
and a, and We horizontal and vertical wind gust velocities. 
Table I lists the stability derivative values used in this 
Simulation. The equations above incorporate the following 
Simplifying assumptions: 

1, The vehicle has been idealized as a rigid airframe to 
which is attached a rotor. 

2. The rotor is described by its tip path plane whose 
orientation determines propulSive and aerodynamic forces 
and moments. 

3. No rotor degree of freedom is considered other than 
control inputs which serve to describe instantaneous 
tip path plane orientation. 

4. Coupling between longitudinal and lateral motion has 
been ignored. 

3. Motion has been linearized as small perturbations about 
a predetermined reference flight path, 

The hybrid computer used consisted of a Scientific Data 
Systems SDS 9300 digital computer, A Comcor CI-5000 analog 
computer, and an Adage AGT-10 stroke written graphics 
processor, as illustrated in Figure 1. The digital computer 
was used to control the analog and graphics processor and 
to generate performance data, A program listing is found in 


Appendix A. In addition to the two display modes to be 
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described, provision has been made for the inclusion of 
flight director symbolism in a third display mode for future 
study. The analog computer gave solutions to the equations 
of motion in real time (Appendix B), and the graphics processor 
displayed textual and graphic information on a cathode ray 
tube in front of the test subject, Pilot control inputs 
were made directly to the analog using electrical voltage 
Signals. The spring restrained cyclic stick operated a gear 
driven potentiometer and was linear in both force and dis- 
placement, The collective lever was unrestrained except for 
a friction brake and also used a gear driven potentiometer 
to generate a signal linear in displacement. The friction 
brake required a slight breakout force. Figure 2 shows the 
physical arrangement of the apparatus. 


Atmospheric turbulence was modeled by a sum of sine 


waves: 
+) 

o = «ama ae where A, = 4,472 Oo, = 0.070 
Ao = 3.536 Wo = 0.210 
A, = 2,236 We = 0.350 
A, = 2,738 wo, = 0.768 
A. = 2.236 Os = 1.765 

and similarly, 

We = E Aysino,t, where A, = 4,472 W4 = 0.420 
Ay = 3.536 gy = 1.260 
Ag = 2.236 Wo = 2.100 
A, = 2.738 as 4,608 
As = 2.236 Oe = 10.050 





The amplitudes and frequencies above were chosen so that 

the integration over w of the power spectral densities 

closely approximated that for the turbulence spectra shown 

in Table II (Ref. 4). The advantage of this turbulence 

representation is that it is repeatable, has a mean value 

of zero, and appears completely random to the test subject. 
_ Two different displays were presented to the test 


subjects for evaluation. The first was basically a real 
‘0 


—— 


world display similar to current instrument installations and 
consisted of a horizon bar and aircraft reference for pitch 
control along with a glideslope deviation indicator. Air- 
speed deviation was provided by angular displacement of the 
aircraft reference "wings." The second display had symbolism 


identical to the baseline display but was quickened with rate 


information according to the equations: 


O guick = 0 + k30 
Douick eae sae ; 


The quickening gains, Ks and ky, were Simple pure gain 


approximations to the pilot transfer functions obtained in 
the modeling procedure. The display, shown in Figure 3, was 
6 1/2 inches wide and 7 inches high, The nominal eye-to- 
display distance was 21/2 feet. 

After considerable informal training, each subject 
received a formal training session in which he used the base- 


tine and quickened display for 90 second runs each. In the 
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formal data sessions, each subject used each display for ten 
90 second runs. The data sessions went as follows: 
The baseline was flown for seven runs, followed by the 
quickened for seven. Then the quickened was flown for 
three runs followed by the baseline for three. Five of 
the "best'' runs for each display were then selected. For 
ten data runs with each display, mean scores were computed 
for each of the root mean square performance measures 
Ge, a) Guaulae. Oy: o)- Then each RMS score for each display 
and run was normalized by dividing it with the respective 
inean RMS score. The normalized RMS scores were then summed 
to provide a single scalar index of performance for each 
display and run. 
The data in Table IV represent the means and standard 
deviations of the RMS scores for the runs with the five 
lowest indices cf performance for each display. 
Figures 5 and 6 graphically portray the more pertinent 
performance data of Table IV for subject A (non-pilot) and 
subject B (pilot), Also indicated are the model predictions 


which best correlate with the data, i.e. F = 0,25. 


IV. OPTIMAL PILOT MODEL AND PREDICTIONS 


The analytic pilot model is shown schematically in 
Figure 4, and is basicaliy the same as that offered by 
Kleinman and Baron (Ref. 5). This optimal pilot model 


includes the following assumptions: 
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1 Subject to inherent limitations, the well trained and 
motivated pilot benaves in an optimal manner. 

2. With certain modifications, the solution to an optimal 
linear control and estimation problem will predict the 
pilot's control characteristics. These modifications 
are: 

a. A pure time delay is included in each of the pilot's 
control outputs. 

b. Each neuromuscular output system may be modeled as 
asians heorder. bag. 

c. Each variable that the pilot observes is assumed to 
contain pilot induced additive white noise which 
scales with the variance of the observed variable. 

d. Each control output is assumed to contain pilot 
induced additive white noise which scales with the 
variance of the control. 

e. The pilot is able to perceive the first derivative 
of a variable, also noise contaminated, but no 
higher derivative. 

The index of performance for the optimization procedure was 

chosen subjectively by the analyst to represent the task 

objectives as perceived by the pilot. Model performance 
predictions are given in Table III for each of the displays. 

The parameter F represents the fraction of the pilot's 

attention devoted to the control task as a whole. For example, 

F = 0.25 suggests that 25% of the vilot's attention is devoted 

to control and 75% is dedicated to non-control activity sucn 


as communication, etc, 
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V. RESULTS 


Table IV lists the root mean square performance data. 
Figures 5 and 6 give a graphical comparison of the most 
pertinent performance measures and the model predictions that 
best correlate with the data (F = 0.25). For each subject, 
the mean value is indicated, as well as plus and minus one 
standard deviation. In each case the pilot model predicted 
performance with a fair degree of accuracy. Additionally, 
the model predicted improvement for each performance measure 
using the quickened display, and the test subjects verified 
this prediction. However, the magnitude of performance 
improvement was greater in each case than was anticipated. 

It is possible that the task as perceived by the pilot may 
not be the same for each display mode. The task as understood 
may be to keep the deviation symbols within a finite distance, 
say 1/8", of their nominal position, and indeed there may 
exist a threshold of error below which the pilot makes no 
control input at all. Since the gains driving the deviation 
symbolism vary considerably between modes, each display then 
defines a new task. One possible arrangement would be to 
provide additional symbolism for the task variables, i.e. 
airspeed and glideslope deviations, which would be common 
between modes. 

The dramatic reduction in cyclic motion, Sys was not 
predicted by the model, Typical time histories of cyclic 
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motion, Figures 7 and 8, show the nature of this improvement. 
Again the task defined by the display mode, in this case 
pitch attitude control, may be different. For both subjects, 
glideslope deviation control was poorest in comparison with 
model predictions, which suggests that more pilot attention 
was given to attitude/airspeed control than to the glideslope. 
Furthermore, collective motion was less for both subjects 
than was predicted. While the cyclic input was of a continuous 
nature, collective inputs were a series of steps (Figures 7 
and 8) in each case. Such control action probably represents 
a wait and see technique used by pilots as a result of the 
lag in collective response. This phenomenon may in part 
explain the discrepancy between predicted and measured 
glideslope deviations since the optimal control model did not 
address this type of control motion. 

Pilot comments indicated that the quickened display 
was more fatiguing than the baseline and represented an 
increase in pilot workload. Unfortunately, no measure of 
pilot workload, such as number of control inputs per minute, 
was recorded. Furthermore, no objective measure of the "weil 
trained and motivated" assumption was available for each 


Subject. 


VI, CONCLUSIONS 


The use of quickening in the cockpit display significantly 
reduced contrcl motion and improved pilot performance. The 


predictive power of the optimal pilot model seemed adequate 
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even though this report represents only a first iteration. 
The quantitative quality of the model should improve as 
adjustments are made. The model index of performance can 

be adjusted to account for the pilot's apparent greater 
concern with pitch attitude/rate deviations as compared with 
altitude/rate deviations. In addition, the displays could 
be revised to eliminate any variation in the task presented 
to the pilot. As the control theoretic pilot model is 
refined and validated through testing it should become a 
valuable tool in the development of cockpit information 


systems and the expanded utility of the helicopter. 
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TABLE I 


UH-1H LONGITUDINAL DERIVATIVES AS USED IN SIMULATION 
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TABLE III 


MODEL PERFORMANCE PREDICTIONS 


Baseline Display Quickened Display 
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TABLE IV 


SIMULATION PERFORMANCE DATA 
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AB 
ACREF 
ANGL 
AS 
CAN 
E16 
E17 
E21 
E23 
E24 
E25 
FRAME 
GS 
GSD 
GSLP 
HORIZ 
IDEV 
IDM 
IER 
IGD 
ITD 
ITEXT 


MODE 


PROGRAM VARIABLES 


Incremental number of display updates during simulation. 
Vector of coordinates defining aircraft reference, 
Angle defining airspeed deviation. 

Airspeed deviation. 

Cosine of given angle. 

Setting for potentiometer O16. 

Setting for potentiometer 017. 

Setting for potentiometer 021. 

Setting for potentiometer 023. 

Setting for potentiometer 024. 

Setting for potentiometer Q205. 

Vector of coordinates defining display frame. 
Glideslope (height) deviation. 

Glideslope deviation rate. 

Vector of coordinates defining glideslope reference. 
Vector of coordinates defining horizon reference. 
Graphics terminal in use, 1 or 2. 

Graphics instruction, draw or move. 

Error type code, 

Graphics instruction, size of graphics input. 
Graphics instruction, size of text input. 

Vector of text input. 

Display mode, 1 or 2. 


Pitch angle. 
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Q Parven Tale. 

RMAS Root mean square airspeed ceviation. 

RMDB Root mean square cyclic movement. 

RMDC Root mean square collective movement. 

RMGS Root mean square glideslope deviation. 

RMTD Root mean square pitch rate, 

RMTE Root mean square pitch angle, 

SAN Sine of given angle, 

SAS Mean square airspeed deviation. 

SDB Mean square cyclic movement. 

SDC Mean square collective movement. 

SGS Mean square glideslope deviation. 

T Time increment. 

THD Mean square pitch rate. 

THE Mean square pitch angle. 

UDR Update rate for the 90 second period. 

UG Vector of horizontal turbulence magnitudes. 
UG1 Specific horizontal turbulence value. 

WG Vector of vertical turbulence magnitudes. 
WG1 Specific vertical turbulence value. 

V Internal clock count, 60 cycles per second. 
X Running variable, abscisséa. 

Y Running variable, ordinate. 

A1,A2,...A9 Sine wave magnitudes. 

OMU1 ,OMU2,..,OMU5 Horizontal turbulence frequencies. 
OMW1 ,OMW2,...OMW5S Vertical turbulence frequencies. 
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Sue ous,,.,oU5 Sine values for horizontal turbdulence. 


SW1,SW2,...SW5 Sine values for vertical turbulence. 
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